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Abstract 


This  paper  will  consIJer  the  nc.ncrical  calcu- 
lation ol  th«  steady  state  rlccatl  equation  util- 
izing a continued  fraction  approach.  Tue  method 
considered  here  has  n.iny  similarities  with  the  par- 
titioned algorithm  approach  [1.2]  considered  by 
D.G.  Lainiotis  for  the  solution  of  riccati  type 
equations.  First  the  matrix  riccati  equation  in 
steady  state  is  converted  (by  a coordinate  trans- 
formation) into  another  riccati  equation  which  can 
be  expanded  In  a matrix  continued  fraction  form.  The 
transformed  matrix  riccati  equation  may  not  be  In- 
vertible .and  therefore  restrictions  are  made  on  the 
general  matrix  case  to  guarantee  the  proper  inver- 
sion properities. 

The  scalar  riccati  equation  is  then  considered  In 
ar.  effort  to  give  incite  into  the  solution  of  the 
matrix  equation.  A numerical  procedure  is  outlined 
to  iterate  on  the.  solution  of  the  scalar  steady  state 
equation.  Upper  and  lower  bound  approximations  arc 
obtained  for  the  steady  state  case  based  on  M terms 
of  the  continued  fraction  expansion.  Theorems  are 
proved  on  the  convergence  properties  of  the  continued 
fraction  expansion  utilizing  mathematical  induction. 

By  using  the  continued  fraction  algorithm,  the  diff- 
erence (or  error)  between  the  riccati  solution  and 
the  first  N terms  of  the  continued  fraction  expansion 
can  be  bounded  from  above  and  below.  ■ 

A theorem  from  N'umber  Theory  demonstrates  that  the 
riccati  equation  is  a quadratic  surd  and  it  must  ex- 
hibit periodicity  in  its  continued  fraction  expan- 
sion. This  result  was  demonstrated  in  [_lj using  the 
partitioned  algorithm  method.  Thu  typo  of  periodicity 
obtained  here  is  slightly  different  but  has  many  of 
the  same  similarities  as  obtained  in  ClJ . 

If  the  problem  of  Interest  were  to  approximate  the 
numerical  solution  of  tb.e  riccati  equation  using  a 
mathematical  appro.ach  called, "Approximation  by  Conver- 
gents",  an  investigation  is  made  on  this  approach  of 
continued  fractions  in  the  calculation  of  an  irrational 
number  ( such  asTl  ) . By  considering  the  continued 


fraction  cxp.aiialon  of  Ti  using  the  l<irgrst  reciprocals 
of  the  difference  bciwc.-'n  Tf  and  its  respective  con- 
vergences as  the  elemi.iil:,  of  the  fraction  expansion, 
this  method  is  most  efficient.  For  the  case  of  the 
rlccatl  equation,  the  periodic  property  of  the  con- 
tinued fraction  coefficients  gives  convergents  which 
are  identical  to  the  periodic  terms  (hence,  most 
efficient).  It  Is  hoped  that  these  results  can  be 
extended  to  the  matrix  case  to  obtain  the  most 
efficient  expansion. 
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THE  MOST  EFFICIENT  METHOD  TO  NUMERICALLY  COMPUTE  THE  SCALAR  SOLUTION  OF  THE 

STEADY  STATE  RICCATI  EQUATION 


DANIELW.  REPPERGER 
Wright  Patterson  Air  Force  Base 
Dayton,  Ohio  (U.S.A.) 


Abstract 


This  paper  will  consider  the  numerical  calcu- 
lation of  th«  steady  state  rlccati  equation  util- 
izing a continued  fraction  approach.  The  method 
considered  here  has  many  similarities  with  the  par- 
titioned algorithm  approach  [l,2]  considered  by 
D.G.  Lainiotis  for  the  solution  of  rlccati  type 
equations.  First  the  matrix  rlccati  equation  in 
steady  state  is  converted  (by  a coordinate  trans- 
formation) into  another  rlccati  equation  which  can 
be  expanded  in  a matrix  continued  fraction  form.  The 
transformed  matrix  riccati  equation  may  not  be  in- 
vertible and  therefore  restrictions  are  made  on  the 
general  matrix  case  to  guarantee  the  proper  inver- 
sion properlties. 

The  scalar  rlccati  equation  is  then  considered  in 
an  effort  to  give  incite  into  the  solution  of  the 
matrix  equation.  A numerical  procedure  is  outlined 
to  iterate  on  the  solution  of  the  scalar  steady  state 
equation.  Upper  and  lower  bound  approximations  are 
obtained  for  the  steady  state  case  based  on  N terras 
of  the  continued  fraction  expansion.  Theorems  are 
proved  on  the  convergence  properties  of  the  continued 
fraction  expansion  utilizing  mathematical  induction. 

By  using  the  continued  fraction  algorithm,  the  diff- 
erence (or  error)  between  the  riccati  solution  and 
the  first  N terms  of  the  continued  fraction  expansion 
can  be  bounded  from  above  and  below. 

A theorem  from  Number  Theory  demonstrates  that  the 
rlccati  equation  is  a quadratic  surd  and  it  must  ex- 
hibit periodicity  in  its  continued  fraction  expan- 
sion. This  result  was  demonstrated  in  [_ljusing  the 
partitioned  algorithm  method.  The  type  of  periodicity 
obtained  here  is  slightly  different  but  has  many  of 
the  same  similarities  as  obtained  in  LJ  • 

If  the  problem  of  interest  were  to  approximate  the 
numerical  solution  of  the  riccati  equation  using  a 
mathematical  approach  called, "Approximation  by  Conver- 
gents”,  an  investigation  is  made  on  this  approach  of 
continued  fractions  in  the  calculation  of  an  irrational 
number  ( such  aslf  ) . By  considering  the  continued 
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fraction  expansion  of  II  using  the  largest  reciprocals 
of  the  difference  between  Tf  and  its  respective  con- 
vergences as  the  elements  of  the  fraction  expansion, 
this  method  is  most  efficient.  For  the  case  of  the 
rlccati  equation,  the  periodic  property  of  the  con- 
tinued fraction  coefficients  gives  convergents  which 
are  identical  to  the  periodic  terms  (hence,  most 
efficient).  It  is  hoped  chat  these  results  can  be 
extended  to  the  matrix  case  to  obtain  the  most 
efficient  expansion. 


1.  Introduction 


The  riccati  equation  has  been  extensively  studied. 
The  approach  used  here  will  consider  a continued 
fraction  method  for  numerical  calculation  of  the 
steady  state  matrix  riccati  equation.  The  continued 
fraction  approach  has  similarities  to  the  partitioned 
algorithm  method  Ql,^  . By  using  the  continued 
fraction  approach,  the  main  advantage  occurs  due  to 
the  fact  that  the  rlccati  equation  can  be  approxima- 
ted by  a finite  number  of  terms  of  the  known  matrices. 
If  a truncation  is  made  after  N terms  of  the  contin- 
ued fraction  expansion,  then  an  upper  and  lower  bound 
for  the  difference  between  the  actual  riccati  solu- 
tion and  the  approximation  can  be  obtained  for  the 
scalar  case.  The  difference  between  the  two  bounds 
can  also  be  calculated  which  gives  an  idea  of  the 
ability  of  this  approach  to  squeeze  the  true  riccati 
solution  between  the  respective  bounds.  Some  examples 
are  worked  to  illustrate  the  advantages  of  this 
approach.  In  order  to  apply  this  approach,  it  is  nec- 
essary for  the  development  of  basic  theorems  before 
the  continued  fraction  method  can  be  utilized. 


2.  Conversion  of  The  Riccati  Equation  Into 
The  Continued  Fraction  Form 


The  problem  of  interest  is  the  solution  of  the 
following  matrix  equation: 

AP  + PA^-PRP  + Q-  0 (1) 

The  unknown  matrix  P is  nxn  and  the  matrices  A,R,  and 
Q are  all  nxn  and  knewn.  It  is  initially  assumed  that  i 

R is  positive  definite;  the  remaining  conditions  | 

t 

I 


! 


J 


(such  as  P and  Q both  synnetric  and  poslttvu  defin- 
ite, A canonical  and  stable,  and  the  controllability 
and  observability  condlllons)  will  not  be  stated 
initially  in  an  effort  to  derive  theorems  based  on 
as  few  assumptions  as  possible.  Theorem  1 illustrates 
the  method  to  convert  the  matrix  rlccati  equation 
(1)  into  the  form  in  which  the  continued  fraction 
approach  can  be  used. 

Theorem  1:  The  solution  to  equation  (1)  with  R 
positive  definite  can  be  specified  as  the  solution 
of  the  following  matrix  equation: 

WW-KW-L-0  il) 

where 

P - W R-1  + R-^  (3) 

Proof: 

Substitute  equation  (3)  into  equation  (1)  yields: 

AWR'l  + AR'^a'''  + MR'^a'^  + R'U'^A^ 

-(WR”^+R-lA'r)R(WR“UR”^A'^)  + Q « 0 
or 

AWR~^  + AR'^a’’’  + WR’^A^  + R“^A^A^  -WWR“^ 

-wr'^a’'’  -r"1a3'rwr-i  - r“1a''''rr-1a'^  + Q - 0 

simplifying  yields: 

ww  + [r-Ia’^r-j^w  + [^-ar-1a^r  - qr3  » 0 

Hence,  if  Che  choice  Is  made  of  Che  form; 


R - -R~1a%  + A 

(4a) 

L » QR  + AR-^^A^R 

(4b) 

then  equation  (2;  will  follow. 

It  is  noted  chat  the  matrix  K may  be  singular 
but  if  R is  symmetric,  Q positive  detlnlce  and  sym- 
metric, then  Che  matrix  L will  he  positive  definite. 

To  conclude  Che  study  on  the  matrix  rlccati  equ- 
ation, the  continued  fraction  approach  ( for  Che 
matrix  case)  will  be  investigated  using  equation  f2) 
in  lieu  of  equation  (1).  The  next  theorem  demonstrates 
the  advantage  of  expressing  equation  (1)  in  Che  form 
of  equation  (2)  for  matrix  continued  fractions. 

Theorem  2 

If  W”^  exists  and  K * 0,  then 

W « K + L 

K+ JL 

K + - 


3.  The  Scalar  Rlccati  Equation 


Consider  the  following  equation: 
2 

aP+Pa-Pr+q“0 

with 

Of  q^  0.  and  a^O 
equations  (4a-b)  now  become: 


L .«  qr  + a 

equation  (2)  now  becomes: 

- (qr  + a^)  - 0 


which  has  a solution  specified  by: 


qr  + a^ 


(6) 


It  is  easily  seen  that  U always  has  one  positive 
root  because:  2 

qr  + a > 0 
P is  specified  by: 


Note  that  under  the  standard  assumptions,  if  U>-  0, 
q>0,  a<^0,  and  r^O,  it  then  follows  Chat  P^  0 

always. 

The  continued  fraction  expansion  for  U is  trivial 
in  this  case.  In  order  to  investigate  the  scalar 
rlccati  equation  in  greater  detail,  convergence 
theorems  will  be  proven  based  on  Euclid's  Algorithm 
[3]  with  tlie  proper  notation  used  in  continued 
fraction  expansions. 


4.  Convergence  Theorems  on  The  Continued  Fraction 
Expression  For  The  Rlccati  Equation 


Consider  the  scalar  continued  fraction  expansion 
of  Che  form: 


+ 1 


32+1 


“3 


- +1 


+ __ 


(7) 


The  continued  fraction  notation  of  equation  (7)  will 
be  specified  by: 


- K + LiK  + LI  K + LI  K+  ]-l  J~^  ]~^  (5) 

if  W”3  exists  and  K-0,  then  W « IL]^^^  which  is 
different  from  standard  matrix  square  root  notation. 
Proof; 

Equation  (2)  can  be  written 

L 

W-K  + LW-1  - K+W 

1/2 

Thus  equation  (5)  follows;  if  K - 0,  then  U • L 
as  indicated  above. 

The  scalar  rlccati  equation  will  now  be  considered 
in  an  effort  to  study  the  numerical  algorithm  presen- 
ted here  and  to  gain  incite  into  the  solution  of  the 
matrix  case. 


V 


<C“o  ’ ®1’  ®2* 


— > 


The  following  definition  is  essential  to  discuss  con- 
vergence properties  of  continued  fraction  expansions. 
Definition : 

Let  a^,  aj^,  a2, be  an  infinite  sequence  of  in- 

tegers, all  positive  except  perhaps  a^.  Define  two 
sequences  of  integers  '{h|A  and^K^  inductively  as 
follows:  ^ ^ 


t'_2  ■ h_i  - 1,  hj  “ «!  hj,!  + h^_2  (Sa) 
K_2  " f.  - 0,  Kj  - aj  + Kj.2  t®**) 
for  i^  0 
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Note: 

Uence; 


Kj-a^K^iK^,  K2>  Kj 

Kn>Vi>^ 

‘'l  ^ - 1 


Theorem  3 

For  any  positive  real  number  x 

<^o*  “1-  ---.  Vl*  ■ 

Proof}  if  N~0,  implies: 


»**H-1  ^ ^-2 

-^N-l  ^ V2 


xh_^  + h_2 

xK_i  + K, 


-2 

which  ii  true  from  the  substitutions  (8a-b). 
For  N"1  the  result  is 


«> 


xhg  + h_j 


«Ko  + *^-l 


which  is  true  by  the  substitutions  (8a-b).  Using 
mathematical  induction  it  is  necessary  to  assume 
that  the  results  hold  for  »®N-1*  ' 

To  calculate  the  partial  sum  for  N terms,  it  is 
necessary  to  use: 

^o.®l*  x^ 

(ajj  + l/x)hj,_^  + hj,_2 
(ajj  + l/x)Kjj_^  + Kjj^2 

- ■"  V2^  ^ S-1  * ^J-1 


<“nVi  ‘S.-2>  S--1 


%-l 


Proof!  ^a^.a^. 


Vr*>' 


*Nj-1  ■*■  ^4-2 


*S-1  ■*■  ^1-2 

from  theorem  (3).  Now  replace  x by  a^^.  This  implies 


'n  ” ^o’®l’  ’ 


“n  "n-i 


+ h, 


'N-2 


“n  ‘Sj-1  *^-2 


(9) 


but  the  definitions  (8a-b)  imply  that  equation  (9) 
can  be  written 


s 


-j,  It  Q.  K.  0. 

The  next  theorem  examines  convergence  of  the 


continued  fraction  expansion  from  the  calculation 
of  an  upper  and  a lower  bound. 

Theorem  5 

As  a lower  bound 

rj  for  J even  forms  an  increasing  sequence,  i.c. 


As  an  upper  bound 

r^  for  J odd  forms  a decreasing  sequence,  l.e. 

ri>r3>r5>r2> 

and  lim  r exists  and  for  every  J ^ 0 
N— • ** 

r^j^lim  rjj  < (j  an  integer) 

l.e. 

^o<'2<"4<'  --  < “Vn<“- V5<^3<'^1 

In  order  to  prove  theorem  5,  the  following  two 
lemmas  must  be  shpwn: 

Lemma  (1) : The  following  equations  hold  for  i^l 


“i^'l-l  - *>1-1  •'i 
which  implies 

•^i  - ’^1-1  ■ 

Lemma  (2)  t for  1*2: 


(-1) 


i-1 


(-1)^-1 
■^i  ■'i-l 


hiKi-a 

which  implies: 


hi_2  Ki 


(-1)*  a. 


(-1)  a^ 


^i-2 


(10a) 


(10b) 


(11a) 


(11b) 


Thus  if  the  results  hold  for  N-1,  then  they  hold  for 
N.  Since  they  hold  for  N*0  and  N=l,  then  they  hold 
for  all  N.  Q.E.D. 

The  next  theorem  is  useful  for  determining  upper 
and  lower  bounds  for  sums  of  the  sequences. 
Definition: 

Define  rjj  = ^ ,aj^, a.j'^  for  all  Integers  N ^0. 

Notice  that  r^  is  the  approximation  after  N terms  of 
the  continued  fraction  expansion. 

Theorcin  4 

r„  - S 


First  to  prove  (lOa-b).  Using  equations  (8a-b) , the 
proof  will  follow  by  mathematical  induction: 

For  i=-l,  it  is  easily  shown  that: 

•»_1  " ’’-2  ’^-1  “ ^ 

Also  using  equations  (8a-b) , the  following  result 
holds  for  1*=0: 

-1 


“-1  Ko 


Now  assume: 


•>1-1  " *>1-2  *^1-1 


(-1) 


i-2 


using  equations  (8a-b)  implies: 

“ *'i-l  *^1  “ ^®i*'i-l'*’*’i-2^*^i-l 

- (hj_j)  (a^Kj_j  + Ki_2) 


®l‘*i-l’^i-l  ■*',*'i-2  '^i-1 


®i’^i-l*’i-l  “ *’i-l  *^1-2 


“ ”^**i-l’^l_2 
since  hj^K 

this  implies 
h 


- hi_2 

i*-!-!  " *'1-1  *^1  ■ 


(-1) 


1-1 


(-1) 


i-1 


1 


(-1) 


1-1 


•^1  *^1-1 


Wl 
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Uencc  Icmna  1 Is  shown.  To  prove  lemaa  C2} i It  is 
observed  that: 


h K , 
o -2 


- h_2  1^0  " *o 


and  equation  (11a)  holds  ior  the  case  !•!  (which  is 
easily  verified) j proceeding  as  before:. 

**1  *^i-2  ” **1-2  *^i  “ ^®i**i-l  * ^I-2^*^l-2 

-hi_2(aiKi_i  Ki_2) 

■ *l''l-l*^l-2  * ’'t-2Ki_2  ~ *t''i-2'^l-l  ■ •*i-2  *^1-2 

“ *1  l*'i-l*^i-2  ” *'1-2  *^l-ll  “ ®i  ^ (-1)  ] 
which  implies* 


“l  ■'i.2  - \-2  "l 


*'<  h , 

— — i-2 


(-1)^  a. 


(-1)^ 

*^1  *1-2 


■^i  " 'i-2 


(-1)'  a. 


*^1  *^1-2 


Hence  the  results  of  lemma  2 hold.  Now  the  proof  of 
theorem  5 will  be  demonstrated: 

Since  each  0,  aj^>0 


(-1)' 

- 


(-D-Ui 

'i~'i-2"  *^i’'-t-2 


WC  fkltv/w  . I J 

(-1) and  (-l)'a^ 

'i  - 'i-i  “ *^1*^1-!  'r'i-2-  ‘^i>'-t-2 

It  follows  for  j“  an  integer: 

r.  ^ r from  equation  (13) . 

2j  \ 2j+2 

From  equations  (12,13)#  we  know  for  aj^^  0 and  j"  an 
Integer: 

'2J-l^'2j-(-l  '2j  <^'2J-1 

Therefore: 

r^<r2<r^<—  and  r^>r^r^ 

It  remains  to  show  that: 

'2n  <^'2n+?<  ^'2n+21-1^21-l 


'2n  <^'2n+2j  ^'2n+2j-1^2J-l 


Tlworem  6 

Tile  continued  fraction  which  represents  a quad- 
ratic surd  is  periodic.  The  term  quadratic  surd  from 
number  theory  (.ij  means  a quadratic  equation  similar 
to  the  riccati  type. 

Proof : 

“ ^o’“l’“2*  •'*n*“o>^l«  »“n*“o'"^ 

Let  - ^o,a2,a2, ,a^ 

Then  x~  implies 

xh  + h , 
n n— 1 

’‘*n  *^-1 

or  X satisfies  an  equation  of  the  form: 

Ax^+¥x  + c • 0 where  A ♦ 0 

because  each  K^^O.  The  converse  of  this  theorem 
can  also  be  shown  but  the  proof  is  more  difficult. 


5.  A Scalar  Example 

Example  1: 

Consider  the  solution  to  the  following  riccati 
equation  by  the  continued  fraction  approach: 

a«  -1,  ■ r«2,  q«l  l.e. 

-p  -p  -2p^  +1  • 0 

2 

or  2p  + 2p  -1  *0  which  has  roots 

1 0.14: 

2-2 


The  positive  root  can  be  written  p«.366 

Now  using  the  continued  fraction  approach,  we  know 
p satisfies  the  following  equation: 

P ■ T - 7 

2 2 

and  w satisfies:  w‘  ■ qr  -f  a 

, -(1)(2)+1  (14) 

w^  - 3 

To  Illustrate  the  applicability  of  this  approach,  if 
w^  - qr  + a^ 
or  w - VT* 

w - 1 + (NJT  - 1)  ^^5) 


w 1 +- 


The  sequence  r^,  T2,  r, — is  monotonically  increas- 
ing and  is  bounded  above  by  r^  and  since  all  r^^O, 
it  has  a limit  In  a similar  manner  the  sequence 

ri,r:,,rc, is  monotonically  decreasing  and  is 

bounded  below  by  r^  and  thus  has  a limit.  The  two 
limits  are  equal  because: 

(-1)^'^ 

lim(  r^  - r ^^1  

i--« 

as  1 — ►-  oc 

Since  each  t*'®  proof  of  theorem  3 is 

complete. 

Before  examples  are  worked  it  is  necessary  to  con- 
sider the  periodicity  of  the  riccati  equation  in  the 
coefficients  of  the  continued  fraction  expansion. 


1 + 1 fT  + V3  -1 1 

2 2 L j 


1 + 1 

2+ 1 

1+ 1 
2+  — 


1 + 


1+  1 

r 1 

2 

' 1 
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Since  NI  - 1.7320508 

Then  the  continued  fraction  of  interest  is: 

1 


1.732 1 + 


1 + 1 


(16) 


2+  1 


Example  2 ^ 

Lets  be  the  irrational  number  U , i.e. 
^ * - 3.14+— 


Hence  a. 


2 + 

Now  it  is  necessary  to  calculate  the  and  h terms. 
Table  I Illustrates  this  calculation.  Now  tRe 
concept  of  "best  possible  approximations"  will  be 
considered  for  the  scalar  riccati  equation  and  for 
any  irrational  number.  TABLE  I 


w 

fl 


3.  To  calculate  a^,  compute 


Ti-  -o 


fl" 

tt-  3. 


.14159 


7.06+ 


hi 

Ki 

h-2  = 0 

K.2  = 1 

fi-Fi  -2 

h.,  = 1 

K.i  = 0 

ho  = 1 

Ko  = 1 

1.0 

3 0 = 1 

h,  = 2 

Ki  = 1 

2.0 

1.0 

ai  =1 

hj  = 5 

<2  =3 

1.667 

-.333 

.6667 

32  =2 

h3  = 7 

K,  = 4 

1.750 

.08333 

-.2500 

33  =1 

h4  = 19 

K4  =11 

1.7272 

-.0227 

.06060606 

34  *2 

hs  = 26 

Ks  = 15 

1.7333 

.0060606 

-.0166666 

35  =1 

ho  = 71 

Ko  = 41 

1.7317 

-.001626 

.00434589 

36  =2 

h7  = 97 

K,  = 56 

1.73214 

.000439 

-.0011905 

3?  =1 

ho  = 265 

Kb  = 153 

1.732026 

-.0001167 

.0003188267 

39  =2 

h,  = 362 

K,  = 209 

1.732057 

.00003127 

-.0000354409 

39  “1 

hio  = 989 

Kio  = 571 

1.732049 

-.0000083795 

.000022893 

ai 0 - 2 

Using 


'"n  * for  n > 0 

''n  - 


•"i  ■ ""i  -1^ 


(-1) 


1-1 


* '■’-2=Tq-iq- 


6.  The  "Best  Possible  Approximation" 


Given  an  Irrational  number  f , the  following  contin^ 
ued  fraction  expansion  is  ''best  possible"  t,5j  for 
numerical  computation  of  ^ 1 
Procedure: 

Using  the  notation  Jwhere  a^  denotes  the 

nearest  Integer  smaller  than  tne  irrational  number 
then  a^  (1”0, ) is  computed  as  follows: 


let 


[?o]  , now  proceeding  as  an  algor Ithmj 

fl”  To-^o  “l  ■ CflJ 


Yo-^o  and  a^  ■ 
inductively  It  follows  that 
1 

fl+1  “ fi~^l 

An  example  will  now  be  worked  with  the  irrational 
number^  to  illustrate  this  approach. 


- [F.] 


Now  compute  a^^  «> 
calculate  fusing  double 


7.  In  this  manner  we 

fusing  double  precision  (29  digits)  on  a 
CDC  6600  computer)  the  following  numerical  results) 


'll 


U- 

ir- 

13  ” 

f5  ■ 

u ■ 

T7 


3.141592653589793238462643383279 

7.0625133059310457697930051531 

15.99659440668569411310599874 

1.0034172310133977856369454934 

292.63459101223866070785852178 

1.5758180949841629954527421461 

1.7366595609113341887024369765 

1.3574791573503535151332712575 


» 2.7973658867611542733590927681 
fg  - 1.2541293985649846777332901666 


’10 


3.9350032135077287460064498245 
“ 1.0695150407541735890218572566 
- 14.385376015764735321479897720 


■14 


This  gives  rise  to  a^“3,  a,"7,  a2”15,  33-!,  a^”292, 
aj-l.a^-l,  aj-1,  ag-e,  ag-1,  a^Q-3,  ajj-l,  a^j" 

which  results  In  the  following  partial  fraction 
expansion: 
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T *■ 


^ 1 

»l  -34i 

7+1 

15+J; 

1+1 

29^+1 

1+1 

1+1 

l+I 


2+j^ 

1+1 

3+1 

1+1 

14+ 


The  first  8 terms  are  calculated  and  expanded  ( the 
computations  are  excessive). 

r « 3 “ 3.0 

r°  = 22/7  - 3.142857142857 

r,  - 333/106  - 3.14150943396226415094 

r,  - 355/113  - 3.14159292035 

rT  - 103993/33102  - 3.141592653011 

r,  - 104348/33215  - 3.14159265392142 

r^  - 208341/66317  - 3.141592653467 

r,  - 312689/99532  - 3.1415926536189 

rg  - 833719/265381  - 3.14159265358107777— 

It  is  noted  that  9 terms  of  the  fraction  expansion 
give  rise  to  accuracy  beyond  9 digits.  The  first  4 
terms  gave  an  accuracy  of  7 decimal  places.  It  can 
be  explicitly  shown  LsJ  that  this  numerical 
approach  is  the  most  efficient  method  to  numerically 
calculate  an  irrational  number  to  any  number  of 
decimal  places.  The  proof  will  not  be  given  here. 

Now  this  approach  will  be  worked  for  the  solution  to 
the  rlccati  equation  C14) . 

Example  3 

Expanding  the  square  root  of  3 to  30  decimal  places 
yields: 

f p > MI  - 1.732050807568877293527446341505 
Now  the  calculation  of  the  a.  and  f will  proceed 
as  for  Tf  . ^ ^ 


a^  - [1.732 I - 1 

" ' “ 1.36+— 

XI  - 1 

Hence  _ 

aj  - CfiH  - 1. 

A summary  of  Che  first  11  terms  yields: 

I - 1.3660254037344363479757026605 

- 2.7320508075688992137673133780 

- 1.3660254037843965712518405766 

- 2.7320508075691909097423020388 
«•*  - 1.3660254037838529560257256617 

- 2.7320508075732481355762416877 
t,  - 1.3660254037762821195840238963 
r J - 2.7320508076297576012776604042 
f “ - 1.3660254036708340246350351294 

i « 2.7320508084168328955065207969 
« 1,3660254022021315334827392008 

Which  results  in  a^^-l,  Sj^-l,  a2“2,  aj“l,  a^“2, 

“5-h  


5 ’ ^6  ’ 7 ’ “8  •‘g-x.a.Q-i. 

In  this  manner,  the  "most  efficient"  method  to 
calculate  the  solution  to  the  rlccati  equation  (14) 
is  as  follows  using  the  continued  fraction  approach: 


1+  J 

2+  J 

1+  1 

1+  __ 

It  is  noted  chat  this  expansion  is  identical  to  the 
one  plv>-n  In  equation  (16).  This  result  occurs  due 
to  the  loiluwing  identity: 

1 + NfT  ^ 


which  results  in: 


. S’ 

I j+2 


for  all  j 


a • 1 ( i odd) 
at^,-  2 (i  odd) 


Since  the  square  root  of  3 is  an  Irrational  number 
(as  the  solution  of  any  quadratic  equation  may  be), 
the  "most  efficient"  method  to  numerically  compute 
the  solution  to  this  rlccati  equation  has  been  demon- 
etxated  by  equation  (16).  This  technique  is  numer- 
ically illustrated  here  in  an  effort  to  develop  an 
algorithm  for  the  matrix  case  and  when  the  solution 
to  equation  (2)  is  not  known  aprlori. 


7.  Conclusions 


A continued  fraction  approach  is  used  to  inves- 
tigate the  solution  to  the  steady  state  rlccati 
equation.  Theorems  are  proven  to  demonstrate  an 
algorithm  which  converges  for  the  steady  state 
scalar  case.  Some  examples  are  worked  to  demonstrate 
Che  "most  efficient"  method  to  compute  Che  solution 
to  the  rlccati  equation. 
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